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Temperature Equilibration Rate with Fermi-Dirac Statistics 
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We calculate analytically the electron-ion temperature equilibration rate in a fully ionized, weakly to 
moderately coupled plasma, using an exact treatment of the Fermi-Dirac electrons. The temperature 
is sufficiently high so that the quantum-mechanical Born approximation to the scattering is valid, ft 
should be emphasized that we do not build a model of the energy exchange mechanism, but rather, 
we perform a systematic first principles calculation of the energy exchange. At the heart of this 
calculation lies the method of dimensional continuation, a technique that we borrow from quantum 
field theory and use in a novel fashion to regulate the kinetic equations in a consistent manner. We 
can then perform a systematic perturbation expansion and thereby obtain a finite first-principles 
result to leading and next-to-leading order. Unlike model building, this systematic calculation 
yields an estimate of its own error and thus prescribes its domain of applicability. The calculational 
error is small for a weakly to moderately coupled plasma, for which our result is nearly exact. It 
should also be emphasized that our calculation becomes unreliable for a strongly coupled plasma, 
where the perturbative expansion that we employ breaks down, and one must then utilize model 
building and computer simulations. Besides providing new and potentially useful results, we use this 
calculation as an opportunity to explain the method of dimensional continuation in a pedagogical 
fashion. Interestingly, in the regime of relevance for many inertial confinement fusion experiments, 
the degeneracy corrections are comparable in size to the subleading quantum correction below the 
Born approximation. For consistency, we therefore present this subleading quantum-to-classical 
transition correction in addition to the degeneracy correction. 

PACS numbers: 52.25.Dg, 05.20.Dd, 25.45.-z, ll.10.Wx 



I. INTRODUCTION AND SUMMARY 

We shall calculate the thermal equilibration rate be- 
tween Fermi-Dirac electrons and Maxwell-Boltzmann 
ions in a hot, fully ionized plasma. We shall do so ex- 
actly to leading and next-to-leading order in the plasma 
number density, and to all orders in the electron fugacity, 
thereby providing an essentially exact result for weakly 
to moderately coupled plasmas. We shall work out this 
problem for two reasons. First, the result is new and is 
needed in some applications. Second, our previous treat- 
ment in Ref. [l[ of the plasma stopping power was per- 
formed in great generality, and the basic idea behind the 
dimensional continuation method, a somewhat subtle an- 
alytic tool that we employ, may have gotten lost in all 
the details. We shall use this opportunity of a simpler 
and specific case to treat the method in a pedagogical 
fashion and to explain it clearly 

Physical systems often contain disparate length or en- 
ergy scales. For example, plasma physics involves hard 
collisions at short distances - ultraviolet physics; and soft 
interactions at large distances entailing collective effects 
- infrared physics. The resulting interplay of short and 
long distances produces the familiar Coulomb logarithm. 
For the electron-ion temperature equilibration rate, and 
for other such processes involving disparate scales, it is 
rather easy to calculate the leading contribution, namely 
the overall factor in front of this logarithm. Although the 
order of magnitude of this leading order term can usu- 
ally be obtained from simple dimensional analysis alone, 
calculating the additional dimensionless factor inside the 



logarithm, the sub-leading term, is quite difficult. 

A new method [3[ employing dimensional continua- 
tion has been introduced to deal with such problems, a 
method that makes the computation of the sub-leading as 
well as the leading contributions tractable. This method 
is based on tested principles of quantum field theory con- 
structed over the last fifty years, and it has been used suc- 
cessfully to calculate well measured phenomena such as 
the Lamb shift [3|, often with much more ease than tra- 
ditional methods. Most recently, the method has been 
exploited in Ref. [l[ by Brown, Preston, and Singleton 
(BPS) to provide an extensive treatment of the charged 
particle stopping power in a plasma, the energy loss per 
unit distance dE/dx of the charged projectile. One of the 
topics treated in BPS was the rate at which elections and 
ions in a spatially homogeneous plasma, starting with dif- 
ferent temperatures, come into thermal equilibrium at a 
common temperature. Here we shall extend this work to 
include the case in which the electron fugacity is suffi- 
ciently large that the electrons must be treated with a 
degenerate Fermi-Dirac distribution. This is the case in 
which Pauli blocking becomes of some importance. 

The degeneracy effects that we treat here come into 
play as the plasma temperature is lowered. We shall cal- 
culate the rate for the general case in which the electrons 
are described by a Fermi-Dirac distribution, with no ap- 
proximations being made to this distribution. That is, 
we shall perform the calculation exactly to all orders in 
the electron fugacity 
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where j3 e = 1/T e is the reciprocal of the electron tem- 
perature and fi e is the electron chemical potential. Note 
that we shall always measure temperature in energy units 
so that fi e does indeed have the correct units of energy. 
We shall assume that the plasma is at most moderately 
coupled, which often implies that the degeneracy correc- 
tions are not large. Nonetheless, we shall work out the 
general case since this is just as easy as treating the case 
of only mild degeneracy, and the general case may prove 
to have some application. 

A plasma is seldom formed in thermal equilibrium; for 
example, a laser preferentially heats the light electrons 
rather than the heavy ions. While a non-equilibrium 
plasma will of course eventually thermalize, it does so 
in several stages. First, the electrons rapidly equilibrate 
among themselves to a common temperature T e . Then, 
somewhat less rapidly, all the various ions equilibrate to 
a common temperature T,. Finally, the electrons and 
ions begin the process of thermal equilibration, with the 
electrons delivering their energy density to the ions at a 
rate 
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We shall calculate the rate coefficient C eI in the fol- 
lowing sections. However, it is useful to examine this 
result now, the rate coefficient in the form (|8.9p . which 
we display below for convenience. This expression for C eI 
is exact to all orders in the electron fugacity, but valid 
only in a temperature regime in which the short-distance 
scattering is described by the quantum-mechanical first 
Born approximation: 
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Here uf is the sum of the squared ion plasma frequen- 
cies ujf, m e is the electron mass, A e is the electron thermal 
wave length, and k 2 is the electron contribution to the 
squared Debye wave number, including electron degen- 
eracy effects. The precise definition of these quantities 
is presented in Section [Til but we should note here that 
we employ rationalized Gaussian units, so that the en- 
ergy of two electrons of charge e a distance r apart is 
given by e 2 /47rr. The structure of the first line that ap- 
pears in our result (|1.3p agrees with the previous result 
of Brysk Q when his Eq. (35) is re-expressed in terms of 
our notation. However, Brysk does not obtain the pre- 
cise result (II. 4[) for In A, but rather only an approximate, 
leading-log evaluation of the usual form ln{6 max /& m i n }. 



Moreover, Brysk 0] also does not obtain the second line 
of our result (|1.3[) . This second line does not contribute 
in the Maxwell-Boltzmann limit of very small fugactiy z e , 
but it does a provide a significant first-order correction 
in z e . 

The previous work of BPS [l[ computed the exact tem- 
perature equilibration rate for a weakly to moderately 
coupled non-degenerate plasma. While this general re- 
sult was rather complicated, the high-temperature limit, 
in which the short-distance Coulomb scattering is given 
by its first Born approximation, has a rather simple form. 
The non-degenerate rate coefficient of BPS can be ob- 
tained from Eq. (|1.3[) by taking the small-fugacity limit 
z e —> 0, which gives 
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where the logarithm In An above is the non-degenerate 
limit of Eq. (fT7i|) . The term 
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in parentheses follows from the well known relation (|2.12p 
between number density and fugacity, and it is just the 
square of the non-degenerate form of the electron's Debye 
wave number. The non-degenerate limit of Eq. (|1.4[) is 
accomplished by the substitution n e — » n e o, and express- 
ing this result in terms of the electron plasma frequency 
provides the from, 
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where we have used the relation k 2 A 2 /27t = fi, 2 w 2 /T 2 . 
The rate coefficient of Eq's. (|1.5[) and (II. 7p is just that 
given previously by Eq. (12.12) of BPS 1], an expression 
that was also quoted in Eq. (3.61) in the introductory 
portion of that work. There is, however, one difference 
between Eq. (fl~7) and the result (12.12) of BPS. Namely, 
the correct term —7 — 1 appearing in Eq. Q1.7[) above 
was incorrectly written as —7 — 2 in Eq's. (12.12) and 
(3.61) of BPS because of a transcription error in passing 
between Eq's. (12.43) and (12.44) of BPS. We have taken 
the opportunity here to correct this mistake. 

Equation (jXTTJ) gives the precise definition of the 
"Coulomb log" (in the non-degenerate limit) for the tem- 
perature equilibration process we have been discussing, 
including the correct additional constant terms, namely 
the terms ln8— 7— 1. We should, however, emphasize two 
points. First of all, the "Coulomb log" is by no means a 
universal quantity, but rather its precise form is process 
dependent. For example, the "Coulomb log" for elec- 
tron conductivity differs from the "Coulomb log" for the 
electron-ion thermal relaxation rate given here Sec- 
ondly we should emphasize that the "Coulomb log" for 
the relaxation rate does not depend upon the ion tem- 
perature Tj but only upon the electron temperature T e . 
Some authors Q incorrectly replace the squared electron 
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Debye wave number k 2 by the fully screened Debye wave 
number k 2 d = k 2 + k 2 . This incorrectly includes the ion 
contribution k 2 , and thus introduces a spurious depen- 
dence on the ion temperature Tj. 

If the plasma temperature becomes very low, then the 
population of bound states must be taken into account. 
In such a regime, our assumption that the plasma is fully 
ionized and weakly coupled breaks down, and our calcu- 
lation is no longer reliable. Hence we shall discuss and 
examine in detail only the limit in which the degeneracy 
corrections are mild, the regime in which the tempera- 
ture is not too low and the electron fugacity is not too 
large. In this limit, only the first-order fugacity contribu- 
tion from the general result (|1.3|) need be retained, and 
therefore the second line in Eq. ()1.3|) . the term omitted 
by Brysk [J], makes an essential contribution. As the 



temperature is lowered, however, the subleading contri- 
bution to the Born approximation becomes comparable 
to the degeneracy correction, and it too must also be ac- 
counted for. In Appendix [Bl we therefore extract this 
subleading correction in the transition region between 
quantum and classical scattering from the general result 
given in BPS. In the following, these two types of cor- 
rections will be called the degeneracy correction and the 
(first) quantum-to-classical transition correction. 

Working to leading order in the fugacity, we only need 
to expand Eq. (|1.3p to linear order in z e — e^ cAle , a result 
contained in Eq. (|8.15p . The quantum-to-classical tran- 
sition correction is contained in Eq. (|B33j) . Expressing 
the fugacity correction in terms of the density according 
to Eq. (|2.18p . the rate coefficient C eI reads 
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where the numerical values of the zeta-function and its 
derivative are 

oo ^ 

((3) = ^^ = 1.20205..., (1.9) 
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and 
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('(3) = — \nk= -0.19812- •• . (1.10) 
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We also write Zi as the ionic charges in units of the 
electron charge e. Before examining expression (|1.8[) 
in detail, we note that the first line is the leading rate 
coefficient calculated in BPS Q, Eq. JOJ above; the 
second line is the first degeneracy correction following 
from Eq. (|1.3[) ; and the third line is the first quantum- 
to-classical transition calculated in Appendix [Bj 

In the last line of Eq. (|1.8|) . the ratio e H /T e describes 
the relative size of the correction, where 

e » = (£) 2 3^ 13 - 6eV ^ 

is the binding energy of the hydrogen atom. For some 
temperature and number density regimes of interest, the 
second and third lines in Eq. (|1.8p become comparable 



in size. Hence, while our main thrust in this paper is 
degeneracy corrections, we must also take into account 
this first quantum-to-classical transition. 

It is conventional to write the Coulomb logarithm 
as ln{6 max /6 min }, where 6 max is a Debye length long- 
distance cutoff, while b m [ n is a short-distance cutoff that, 
depending upon the circumstances, is either a classical 
distance of closest approach 6 c i ~ e 2 /T e or a quantum 
wave length b qm ~ A e . Often, an interpolation is made @ 
between these two limits by writing 
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Such an interpolation gives a first correction proportional 
to the proper quantum expansion parameter e H /T e , but 
it fails entirely to produce the proper logarithmic be- 
havior (e„ /T e ) ln{e H /T e } displayed in the last line of 
Eq. (EBJ. 

Figures[T]and[2]illustrate the size of the degeneracy and 
the first quantum-to-classical transition corrections in an 
cquimolar deuterium-tritium plasma, for electron num- 
ber densities of n e — 10 25 cm -3 and n e — 10 26 cm -3 , re- 
spectively. The solid curves denote the size of the degen- 
eracy corrections relative to the leading BPS term, the ra- 
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FIG. 1: The relative size of the degeneracy correction and 
the first classical-to-quantum correction as a function of tem- 
perature in keV. The plasma is equimolar deuterium-tritium 
at an electron number density n e = 10 25 cm -3 . These cor- 
rections are relative to the leading non-degenerate BPS rate: 
the degeneracy correction (the solid line) is the ratio of the 
second to the first line in Eq. (|1.3[1 . while the classical-to- 
quantum correction (the dashed line) is the ratio of the third 
to the first. The electron temperature runs between values 
0.1 keV and lOkeV. Our calculation ceases to be valid at low 
temperatures, and this is indicated by the vertical dotted line. 



tio of the second to the first line of Eq. ()1.8|) . The dashed 
curves show the corresponding size of the first quantum- 
to-classical transition, the ratio of the third to the first 
line of Eq. (|1.8|) . Note that both corrections are compa- 
rable in magnitude between these two electron number 
densities. We have plotted the corrections for tempera- 
tures between T e — 0.1 keV and T e = lOkeV. However, 
below about 0.5 keV, ionization and other strongly cou- 
pled plasma effects become important, and our formalism 
breaks down. Also, at these lower temperatures, higher 
order fugacity terms in Eq. (| 1 .4|) become important, and 
this could change the low-temperature behavior of the 
solid curves; therefore, one should trust Eq. (|1.8p only to 
the right of the vertical dotted line at T e = 0.5 keV. For 
densities below n e = 10 25 cm -3 , the degeneracy correc- 
tion is much smaller than the quantum-to-classical cor- 
rection. The situation is reversed for densities greater 
than n e = 10 26 cm -3 , where degeneracy effects dominate 
over the quantum-to-classical corrections. Interestingly, 
the density and temperature range in which the degener- 
acy and the quantum-to-classical corrections are compa- 
rable lies in the regime relevant for inertial confinement 
fusion. 

We turn now to present the details of our work. After 
establishing the conventions and notation that we em- 



FIG. 2: Same as Fig. [Tj except the number density is 
n e = 10 26 cm -3 . Again, the solid line is the degeneracy cor- 
rection and the dashed line is the first quantum-to-classical 
correction. 



ploy, we review and explain in some detail the dimen- 
sional continuation method that forms the basis of our 
calculation, and we relate this method to previous work 
that utilized convergent kinetic equations. After estab- 
lishing this foundation, we then perform the calculation 
of the degeneracy corrections to the rate at which the 
electrons and ions come into thermal equilibrium. Ap- 
pendix [A] contains the details of the plasma dielectric 
function needed in the text, and Appendix |B] extracts 
the first correction in the transition to classical scatter- 
ing from a general formula given in BPS [l[. 

But before passing to these details, we would like to 
make a comment. There are a large number of very able 
theoretical physicists who work on formal developments 
that have no contact with experiment. In the chance that 
one of them may encounter this paper, we would like to 
encourage them to once in a while apply their skills to 
topics such as that treated here. The field would profit 
from their work. In this regard, we quote a penciled poem 
written in 1957 on note-paper embossed with the header 
"The White House, Washington" : 

A fact without a theory is like a ship without a sail, 
is like a boat without a rudder, 
is like a kite without a tail. 

A fact without a theory is like an inconsistent act. 
But if there is one thing worse in this confusing 

universe, 
it is a theory without a fact. 
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II. 



CONVENTIONS AND NOTATION 



Changing to dimensionless variables, 



We will treat the ions with Maxwell-Boltzmann statis- 
tics and the electrons with Fermi-Dirac statistics. The 
thermal equilibrium form of an ion phase space density 
fi(pi) thus reads 

/i (p i ) =e -ft(s i (p,)- W ) ) (2.1) 

while for electrons in thermal equilibrium, 

1 
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Letting the index b refer to either the ions i or the elec- 
trons e, the kinetic energy is 



E b {p b ) 



(2.3) 



and the inverse temperature and chemical potential are 

b = l/T b and ^ • (2.4) 

Since we shall work in an arbitrary number of dimensions 
v, each species number density appears as 
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where $b is spin-degeneracy factor. For electrons g e = 2. 
We are using a notation for the distribution functions 
fb in which the species index b implicitly includes spin 
degrees of freedom. It is inconvenient, however, to use 
this notation for the number density rib- Measurements of 
the species density are usually insensitive to spin degrees 
of freedom, and we shall therefore denote the number 
density of the species (including all the spins) by rib- This 
accounts for the factor of Qb in Eq. (|2.5|) . For the ions 
in thermal equilibrium, the integral (|2.5p is a product of 
trivial Gaussian integrals, and so 
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where we define the thermal wavelength for species b as 

1/2 
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For the electrons in thermal equilibrium, we first pass to 
hyper-spherical coordinates to write 
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Here £l„_i is the area of a unit hyper-sphere in a space 
of v dimensions that is evaluated in the next section with 
the result 
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allows us to express the electron number density as 
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When the quantity —f3 e fi e becomes very large (and 
positive), Fermi-Dirac statistics pass to the Maxwell- 
Boltzmann limit, with the denominator above becoming 
a simple exponential. In this limit, the x-integration be- 
comes the standard representation of T(v/2), and we see 
that the number density in this small fugacity limit is 
given by the Maxwell-Boltzmann form (|2.6[) . as it must 
be. Expanding Eq. (|2.11[) to second order in the fugacity 
z e = e 13 "^" and using Q e — 2 gives, for the the physical 
case of three dimensions, 
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Note that the first correction, which decreases the num- 
ber density, is simply the fugacity divided by a numeri- 
cal factor of order unity. By small fugacity, we therefore 
mean that e^^' < 2 3 / 2 ~ 2.8. 

As we shall see in the following section, in a space of v 
dimensions, the energy of two charges e a distance r apart 
is proportional to e 2 /r u ~ 2 . Since the units of a number 
density n are (length)"", we conclude that e 2 n has the 
units of energy over length squared, independently of the 
spatial dimension v. In particular, 



e 2 n b 



(2.13) 



is the squared plasma frequency for species b with the 
fixed dimension of an inverse-time-squared, regardless 
of the spatial dimensionality v. The situation for the 
squared Debye wave number is essentially the same, ex- 
cept that, as noted in Appendix lAl in general, this quan- 
tity is defined in terms of the fluctuations in the number 
density, and so 
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For Maxwell-Boltzmann statistics, the derivative that ap- 
pears here just reproduces the particle density in accord 
with the fact that classical particles are described by 
Poisson statistics. However, for Fermi-Dirac statistics, 
one must use 
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d(fibHb) r(z//2) J x [ e -PbMb e * + l] 
The inequality that appears here implies that 
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For the dilute case in three dimensions, including the 
first correction in the fugacity, the electron Debye wave 
number is given by 



v = 3 : 
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In the first-order fugacity correction that appears here 
we can use the lowest-order result 



(2.18) 



to compute the fugacity and thus write 
v = 3 : 



Pe e n e 
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III. METHOD 

Since the method of dimensional continuation that we 
shall use is novel and perhaps subtle, we present here a 
pedagogical account of its basis. 



A. Disparate Length Scales; Expansion Parameter 

The electron-ion energy exchange brought about by 
their collisions in a plasma involves a Coulomb interac- 
tion that is Debye screened at large distances and, as we 
shall see in the course of our work, cut off at short dis- 
tances by quantum effects. As we shall sketch below in 
Subsection IIII CI for arbitrary spatial dimensions v, the 
familiar elementary description of this energy transfer for 
v = 3 dimensions involves the impact parameter integral 



db f 6 ma>: 
y =ln 



(3.1) 



Here 6 m i n is the minimum distance of closest approach 
that, in the quantum limit that is relevant here, is set 
by the scale of the electron thermal wavelength A e . That 
is, &min is some numerical constant times A e . The up- 
per limit on the impact parameter integral is set by the 
electron Debye length k~ , with 6 max some numerical 
multiple of k" 1 . Thus 



(3.2) 



The purpose of the dimensional continuation method is 
to precisely determine the numerical constants that ap- 
pear here. 

Our method applies when the ratio b max /b m - ln is large: 
In this case, the dimensionless parameter X e K e is small, 



and we shall use it as our expansion parameter. As we 
noted in the previous section, the Debye wave-number 
K e always has the dimensions of an inverse length, even 
at arbitrary spatial dimension v. Hence, X e K e is a con- 
venient parameter to employ in our dimensional contin- 
uation scheme because it remains dimensionless as the 
number of spatial dimensions v is varied. Moreover, as we 
shall see, it is the combination X e K e that directly arises 
as our computations progress. 

The electron plasma coupling strength is characterized 
by the ratio of the Coulomb electrostatic energy of two 
electrons a Debye distance apart divided by the temper- 
ature. In the physical space of three dimensions, v = 3, 
this is the dimensionless parameter 
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The perturbative expansion of plasma thermodynamic 
parameters involve a series of ascending integer powers 
(up to additional logarithmic corrections) of the coupling 
constant g e . Except for different conventions that can al- 
ter a trivial overall factor, the electron quantum Coulomb 
parameter is defined as the Coulomb energy for two elec- 
trons a thermal wave length apart divided by the tem- 
perature. For three dimensions, this reads 



v = 3 : rj e = {3 e 



Hence, in three dimensions, 



(3.4) 



(3.5) 



Thus our expansion parameter \ e n e is essentially the 
plasma coupling parameter g e , albeit divided by the 
quantum parameter rj e . Accordingly, one could equiv- 
alently work in terms of the coupling g e as we have done 
in the past [l|, but here it is more convenient to use A e K e , 
and so this we shall do. 

Our work applies to fully ionized plasmas where the 
temperature is large and thus the parameter rj e is small. 
The condition that \ e n e be small requires that the 
plasma coupling g e be even smaller than n e . To put 
this in perspective, we recall that even if Fermi-Dirac 
statistics are required, the Debye wave-number is smaller 
than that given by the Maxwell-Boltzmann form with the 
same temperature and density. Hence we have 



a: 



2 kI 



< KP, 



e 2 n e , 



(3.6) 



where the electron number density on the right-hand side 
of this equation is given by Maxwell-Boltzmann statis- 
tics. Using now the number density (|2.6p in the Maxwell- 
Boltzmann limit and the definition (|2.7p of the thermal 
wavelength, we find that 
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Thus, even for somewhat large electron fugacities e /3 " tJ ' e , 
the expansion parameter X e n e will be small provided the 
temperature is reasonably large. 

To explain further the utility of A e the appropri- 
ate expansion parameter, we examine the situation when 
the spatial dimension v departs from its physical value 
v = 3. In this case, as we shall soon see in the next sub- 
section, the Coulomb potential a distance r away from 
a point charge has the dependence r~("~ 2 '. Thus the 
plasma coupling and quantum Coulomb parameters have 
the form 



a 2 v—2 



and so g e /rj e ~ (A e K e )" 2 or 



A e K e ~ | 



l/(v-2) 



(3.8) 



(3.9) 



This emphasizes that although the form of the coupling 
A e K e that we employ here does not change as the spa- 
tial dimension is altered, its form in terms of g e /Ve does 
depend upon this dimensionality. 



B. Idea of Dimensional Continuation 

We have already seen explicitly how a geometrical 
quantity, namely the number density, can be computed 
in a space of arbitrary dimensionality v. In fact, all fun- 
damental theories can be formulated in a world that has 
space of arbitrary dimensionality. Modern quantum field 
theory, the mother of all physical theory, is generally for- 
mulated for spaces of arbitrary dimensionality in order 
to regulate it. [See, for example, Ref. [8j.] The well 
known BBGKY hierarchy of coupled equations that de- 
picts general kinetics can obviously be written in a space 
of arbitrary dimensionality v. For v > 3, the Coulomb 
force acts as a short-range force; for v < 3, it acts a long- 
range force. Although the complete BBGKY set of cou- 
pled equations for Coulomb forces must remain valid for 
arbitrary spatial dimensionality it cannot be approx- 
imated by the Boltzmann or Lenard-Balescu equations 
for general v values. To leading order in the density, the 
Boltzmann equation describes the short-distance, hard 
scattering correctly while the Lenard-Balescu equation 
correctly describes the long-distance, collective interac- 
tions. Hence, to this leading order, the BBGKY hierar- 
chy of equations reduces to the Boltzmann equation for 
v > 3, but for v < 3, the BBGKY hierarchy reduces 
to the Lenard-Balescu equation. We shall see how this 
works out in detail as our work progresses. 

Here we introduce the idea of dimensional continua- 
tion by examining the case of electrostatics. The Poisson 
equation for a point charge Ze in v dimensions reads 



V 2 (l/) (r) = Ze8 {u) {v). 



(3.10) 



Its solution may be expressed as a Fourier integral: 

d v k Ze 



^)(r) = 



(2n) v k 2 



(3.11) 



As it stands, this integral is defined for all positive inte- 
ger dimensions v. As is well known from the theory of 
complex functions, an analytic function is defined from 
its values on the positive real integers, provided that the 
function does not diverge rapidly at infinity in the com- 
plex plane [2j. But for our equations, we can obtain this 
extension by explicit calculations. For the case at hand, 
we first write 



1 

k 2 



1 



ds t 



(3.12) 



and interchange integrals to encounter 
d v k 



(27T)" 

d v k 



„ — k s ik-r 

e e 



(27T)" 



exp 



k ?:— 

2s 



exp< — — 



4ns 



v/2 



(3.13) 



Here we have completed the square and used the variable 
in the square brackets as the new integration variable to 
obtain a product of v ordinary Gaussian integrals whose 
evaluation produces the final result. The change of vari- 
ables from s to t = r 2 /4s now gives 



(r) = 



Ze 



Ze 

-v-2 



v/2 



v/2 







dt 

T 



U t t {v-2)/2 e - 



If 

4 



u-2 



(3.14) 



since the t integral is a standard representation of the 
gamma function. This result now defines an electrostatic 
potential for any value of v in the entire complex plane. 

As a mathematical application of this result, we note 
that it gives the electric field 



E(r) 



Zev 



v/2 



(3.15) 



Hence Gauss' law applied to a sphere of radius r, 

S{r) r • E(r) = Ze , (3.16) 

informs us that this sphere has an area S = r v ~ x 
where 



2W 2 



(3.17) 



is the area of a [y — 1) -dimensional unit hypersphere em- 
bedded in the j/-dimensional space. 



8 




FIG. 3: Short-distance or ultraviolet (UV) physics dominates 
in dimensions v > 3. Long-distance or infrared (IR) physics 
dominates when v < 3. UV and IR physics are equally im- 
portant in v = 3. 



There are physical implications that follow from (|3. 14|) 
of the electrostatic potential of a point charge in v di- 
mensions. These are brought out in Fig. [3] As the figure 
shows, the Coulomb potential of a point charge becomes 
more singular at the origin as the spatial dimension v 
increases: the physics at short distances is increasingly 
emphasized as the spatial dimension v increases. Since 
short distances correspond to high wave numbers, this 
is equivalent to stating that large v emphasizes ultra- 
violet physical processes. Conversely, as the spatial di- 
mension v decreases, the potential falls off less rapidly 
at large distances: the physics at large distances be- 
comes ever more important as the spatial dimension v 
decreases. Since long distances correspond to low wave 
numbers, this is equivalent to stating that small v em- 
phasizes infrared physical processes. As we shall see, the 
electron-ion energy exchange can be computed with the 
Boltzmann equation for v > 3 since it correctly accounts 
for hard scattering. The result, however, has a simple 
pole that diverges as v — > 3 from above. Conversely, 
the electron-ion energy exchange can be computed from 
the Lenard-Balescu equation for v < 3 since it correctly 
accounts for the long-range screening brought about by 
the collective, dielectric effects in the plasma. The re- 
sult, however, has a simple pole and diverges as v — » 3 
from below. These general features are brought out in 
the simple computation of the next subsection. 



C. Energy Loss Structure in v Spatial Dimensions 

To illustrate the remarks that we have been making, 
we consider the lowest-order energy loss of an electron 



passing by a fixed point of charge Ze. In zeroth order, 
an electron with impact parameter b and perpendicular 
velocity v, so that b • v = 0, simply follows the straight 
line b + v t as a function of the time t. In first approxi- 
mation, with the electric field E given by Eq. (|3 . 1 5[) . the 
electron acquires a momentum transfer 



Ap = -e 



eftE(b + vi) 



(3.18) 



in passing by the fixed point charge Ze. This gives an 
energy change 



AE 



2m e 



(3.19) 



It is a straightforward matter to check that this gives, up 
to numerical factors of no importance, 



AE 

TO. 



1 fZe 2s 2 



1 



,v-1 



(3.20) 



In v spatial dimensions, an element of cross section is 
given by da = f2„_2 b v ~ 2 db. Hence, the weighted energy 
loss has the form 



da AE 



db 



(3.21) 



This example explicitly demonstrates that large v is 
dominated by short-distance physics and small v is dom- 
inated by long-distance physics. Moreover, it shows ex- 
plicitly that v = 3 is the dividing line between these two 
regions. To bring this out, all we need do is to note that 
for v > 3 the impact parameter integral is not sensitive 
to the large distance cut off, and we may take the limit 



oo to obtain 



v > 3 



I>{v) 



db 
6^2 



Conversely, for v < 3, we may set & m i n = 0, with 



v < 3 



I<{v) 



b 3 - 



db 



b v - 2 Z-v 



(3.22) 



(3.23) 



The results displayed are the dominant forms in the two 
different regions of spatial dimensionality v. 



D. Implementation of Dimensional Continuation 

The situation that we have just described leads to a 
well defined result because it is akin to the following ex- 
ample. Suppose that we have a theory that is well defined 
in the neighborhood of the physical dimension v = 3, and 
that the theory contains a small parameter e. Moreover, 
suppose that we need to evaluate a function F that de- 
pends upon this small parameter e in the following fash- 
ion. For v > 3 the leading behavior of F goes like e 3_l/ , 



9 



analytically continue 
around the v = 3 pole 



F<{v;e) 




F<(u;e) 



-(3-*) 



-\u-3\ 



LO: large when e < 1 



NLO: small when e < 1 



FIG. 4: The analytic continuation of F < {v; e) from v < 3 to 
the region v > 3: the same expression can be used for F < (v\ e) 
throughout the complex plane since the pole at v = 3 can 
easily be avoided. Note that the quantity F < (v;e) ~ e 
is leading order in e for i/ < 3. However, upon analytically 
continuing to v > 3 we find that i ?< (^;e) ~ e ' ^ 3 ' which is 
next-to-leading order in e relative to F > (u;e) ~ e~'"~ 3 '. 



(*-3) 



and the function f has a simple pole in ^ as v — ► 3 from 
above. Conversely, for < 3, the leading behavior of 
the function F goes like t v ^ z and the function F has a 
simple pole in v as ^ — > 3 from below. That is, we have 
the leading terms 



v > 3 



F > (z/;e) = A>(i/) e 



3-1/ 



(3.24) 



and 



v<'i: F < {v\e) = A < {v)e v ~ i . (3.25) 

Since the two contributions each have poles in v. 
R > 



and 



A>{v) 



A<{v) = 



+ r> +0(u-3), (3.26) 



+ r K + 0(3 - //) . (3.27) 



The function F < (iv; e) is of leading order in the expansion 
parameter e for v < 3. Since it is an analytic function 
of z/, it may be continued throughout the complex v- 
plane. When it is analytically continued to v > 3 it 
becomes of subleading order. This behavior is depicted 
in Fig. |U Exactly the converse situation applies to the 
function F > [y 1 e). 

Therefore, in the neighborhood of v = 3, 



F(u;e) = F>(v;e)+F<(v;e) 

= A > {v)<?- V + A < (u)e u - 3 



(3.28) 



is accurate to leading and sub-leading order in e. For 
v > 3, the term with the coefficient A > (z/) is dominant 
while that with the coefficient A K {v) is sub-dominant. 
For v < 3 the roles of the the dominant and sub-dominant 
terms are interchanged. It should be emphasized that the 
addition of the two terms contains no double counting 



since in each region one term, and one term only, dom- 
inates. Since the theory is well defined at the physical 
dimension v = 3, the poles must cancel, which requires 
that 



R> = R 

Using 

e^ v " = e- 
we now have, in the neighborhood of v = 3, 



,±(w-3) lnt 



(3.29) 
(3.30) 



F(u;e) = 



D -(i/-3) ln< 



v-3 
2i?> lne 



(3.31) 



We must emphasize that this method of dimensional con- 
tinuation provides not only the coefficient 2R out in front 
of lne (which is often not too difficult to compute), but 
the constant r > + r < in addition to this logarithm (which 
is often difficult to compute). 

A relevant example is provided by the simple model 
of the energy loss presented in the previous subsection. 
According to our general method, in the neighborhood 
of v = 3 we must have 



I{y) = F>{v)+I < {v) 



l3-i/ 



v — 3 3 — v 



(3.32) 



The v — > 3 limit produces 



l(y) 



-In 



(3.33) 



This is precisely the value (|3.ip of the familiar impact 
integral evaluated directly in three dimensions 

Another instructive example of this method is provided 
by an examination of the modified Bessel function K v (z) 
for small v and small z. This is discussed in Ref's. [|[ and 
and in more detail in Ref. 0. These works should 
be consulted if the explanation already given is not con- 
vincing. 

Although we have sketched the basic idea of our dimen- 
sional continuation scheme leading to the result (|3.31|) . 
in this paper we shall apply it in a slightly different form, 
a form similar to the more physical example that led to 
the result (|3 . 33[1 . As we have explained, the electrostatic 
potential in v spatial dimensions has the functional form 



e/i 



so that the energy between two point charges a 



distance r apart is proportional to e 2 /r !y ~ 2 . As we shall 
see explicitly in our work below, the electron-ion energy 
exchange rate contains an over all dimension-bearing fac- 
tor j3 e e 2 . This factor has the dimensions of length to the 
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power v — 2. To define a quantity whose physical dimen- 
sion does not vary as the spatial dimension varies, the 
factor j3 e e 2 must be accompanied by a factor of length 
raised to the power 3 — v, which gives a result that has 
a constant factor of 1/length in all spatial dimensions 
v. For the v > 3 contribution, a scattering term with 
a length cutoff given by the electron thermal wavelength 
A e , the needed dimensional factor is given by A 3 ^ as 
we shall explicitly find below. For the v < 3 contribu- 
tion, a long-distance Debye screened interaction term, 
the needed dimensional factor is given by (1/K e ) 3_l/ as 
we shall also explicitly see below. Thus, in all spatial 
dimensions near v — 3, the rate has the structure 



G(y) = /3 e e 2 



\ A - v B>(y) 



1 



3-k 



B<{v) 



and for v near v = 3. 



and 



B>{v) 



B<{v) = 



R 



^-3 



b > 



(3.34) 



(3.35) 



I? 



Z-v 

Writing B = b > + 6 <5 we find that for v near v = 3, 

R 



(3.36) 



G{v) = (3 e e'X 



P e e 2 



2 \3-i/ 



{i-(A e Ke y- 3 } + 



B 



Rln{\ e n e } + B 



(3.37) 



in which the final line gives the v = 3 limit. 

An objection could be raised that we have not shown 
explicitly that larger sub-leading terms are not present. 
We have extracted terms that have the generic behavior 
e 3 ~ l/ for v > 3 and t v ~ z for v < 3. One might ask if there 
are additional terms with a power law dependence be- 
tween e 3_l/ and e^ -3 . However, simple dimensional anal- 
ysis shows that such terms of intermediate order cannot 
appear. The point is that the physics involves only two 
different mechanisms that dominate at large and small 
scales. These two different mechanisms involve differ- 
ent combinations of basic physical parameters and hence 
give quite different dependencies on the small parameter 
when the spatial dimension v departs form v = 3. 



Liboff, in his Eq. (2.75), writes the transport equation 
for a homogeneous system such as we consider as 



~dt~ BQ 



Lq — R, 



(4.1) 



where Bo is the Boltzmann collision integral, Lo the 
Lenard-Balescu integral, and R is a renormalization term 
that cancels the singularities in Bq and Lo [14J. We 
should note that starting off with admittedlyinfinite, and 
therefore undefined quantities, as in Ref. 13], is at best 
a heuristic procedure. This is to be contrasted with the 
renormalization procedure performed in modern quan- 
tum field theories where the starting point is a rigorously 
defined, finite theory because the starting point is a reg- 
ularized theory. At any rate, the infinite renormalization 
term R is expressed formally as a double integral over 
both impact parameters and Fourier wave numbers. The 
integral over impact parameters b is broken up into a 
large impact parameter part b > bo and a small impact 
parameter part b < bo, R — R(> b ) + Ro{< b ). It is 
then shown that R has a formal construction such that 
both Bq — R(> bo) and Lo — R(< bo) are finite. 

The transport equation is thus rendered finite. Li- 
boff concludes, "So we find that the combination of col- 
lision integrals gives a reasonable model for a conver- 
gent plasma kinetic equation." Our goal, however, is not 
just to find a "reasonable model," but to calculate the 
Coulomb logarithm in a precise and rigorous fashion. To 
leading order in the plasma density, we shall not only 
compute the coefficient out in front of the logarithm, but 
also the constants that appear in addition to the loga- 
rithm. See Ref. [2j for more details. 

Although Gould and DeWitt [ll| also separate the 
right-hand side of the transport equation into three 
terms, they do so in such a fashion that each of the terms 
is finite and well defined. As shown in Appendix B of 
BPS [l[, the formulation of Gould and DeWitt correctly 
gives the constant term as well as the leading Coulomb 
logarithm, and as far as these terms are concerned, their 
work is mathematically equivalent to our method of di- 
mensional continuation. Both are accurate to 0(g 2 ) in 
the plasma coupling, and no better. The trouble with 
their formulation is that it also produces a subset of 
higher order terms, and there is no reason that these 
provide a more accurate evaluation. As is well known, 
the inclusion of partial subsets of higher-order terms can 
sometimes give less rather than more accurate results. 



V. ENERGY AND TEMPERATURE RATES 



IV. CONVERGENT KINETIC EQUATIONS 

A number of authors , [n], E3] have proposed vari- 
ous versions of plasma kinetic equations that have neither 
short nor long range divergences. This work is summa- 
rized in the book of Liboff 13], which we shall outline 
here and then relate to our method of dimensional regu- 
larization. 



The rate of change in the electron energy density trans- 
ported to all the ions species vanishes when the two sub- 
systems have the same temperature. Hence we may write 

^L = - CeI (T e - Tl ). (5.1) 

Since energy flows from the electrons to the ions when the 
electrons are hotter than the ions, C eI is positive. Since 
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the total energy is conserved, the rate at which energy is 
transferred from the ions to the electrons, d£ Ie /dt, has 
the same coefficient C eI but an overall sign change or, 
equivalently, 



d£ Ie 
~dt 



C eI (T, — T e 



(5.2) 



A change in the energy of a subsystem in the plasma 
produces a corresponding change in the temperature of 
that subsystem. Thus, for the electrons, 



while for the ions 



A£ Ie 



E, 



AT,, 



A£ t = Cl AT, . 



(5.3) 



(5.4) 



Here, since the plasma interactions do not change particle 
number densities, the specific heats c e and c 7 are those 
at constant volume. Since A£ eI is an energy density, 
these are the specific heats per unit volume. For a hot 
plasma that is not strongly coupled, the case treated in 
this paper, these specific heats are given by the familiar 
ideal gas results: 



3n e /2 



and 



3n,/2 ; 



(5.5) 



where n, is the total ionic density, the number of all the 
ions per unit volume. Thus Eq's. (|5.ip and (|5.2p are 
equivalent to 



dt 



-7ei (T e - T T ) 



(5.6) 



with j eI = C eI /c e ; and 

dTj 
~~dt ~ 



(5.7) 



with 7 Je = C eJ /cj. Moreover, the rate at which the sep- 
arate temperatures approach one another is given by 



rf(T e -T) 
dt 



-r (r e - t) , 



in which 



(5.8) 



(5.9) 



r = c, ,[- + - 



We turn now to compute the rate coefficient C e7 . 



VI. BOLTZMANN EQUATION: 
SHORT-DISTANCE PHYSICS 

We first work in v > 3 dimensions where the short- 
distance physics dominates. Thus the rate of change of 
the electron distribution is described by the Boltzmann 
equation with Fermi-Dirac statistics for the electrons and 
Maxwell-Boltzmann statistics for the heavy ions. The 
Boltzmann equation for the electron distribution, includ- 
ing the Pauli blocking of the scattered electrons and two- 
body quantum scattering effects, reads 



dfe(Pe) 

dt 



d V V 'e d V p' % d v p % 



T 



J (2txK) v (2nh)» (2nh) 

V / \ ^f/fcg —II// —III, £*l 1 1/ 2 

/e(Pe)/i(pD fl - /e(Pe)l - fe(Pe)fi(Pi) k ~ fe(p' e ) 



(6.1) 



where T is the amplitude for the two-body scattering collision ei —> e' i', and we have omitted the spatial convection 
term on the left-hand side because we are concerned with spatially uniform plasmas. Since the electrons are in thermal 
equilibrium with each other, there is no electron-electron interaction contribution to this time derivative. The electron 
kinetic energy density — the electron energy per unit volume — has the same form as the number density (|2.5p save 
that an additional factor of E e (p) = p 2 /2m e appears in the integrand. Hence the rate at which this energy density 
changes because of the electron ion interactions is given by 



d£g 
dt 



= 2 



d u Pe pi 9/ e (Pe) 



{2-Khy 2m e dt 



(6.2) 



where the factor of 2 multiplying the integral accounts for the electron spin degeneracy. Using the crossing symmetry 
Pe <-> p' e and Pi <-> pj of the scattering amplitude T in (|6.1|) , the rate of energy exchange from the electrons to the 
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ions (|6.2[) can be written as 



dt 



d" P ' e d v p\ d" Pe d v Pi 



J {2irhy (2nhy (2ttK) u (2tt/i) 1 
{2*hy 6"(p' e + p', - p e - pi) (2tt/*) J ( P ^~ Pg 



,/2 



-.'2 



2m; 



2m P 



/e(p e )/i(Pi) l-/e(Pe) 



(6.3) 



where the factor of two in front of the sum is the spin-degeneracy g e — 2 for electrons. We have placed a "greater- 
than" superscript on the left-hand side of the equation since we are now computing the v > 3 contribution. We start 
by performing the p/-integration in Eq. (|6.3p . using the momentum conserving delta- function to set 



P* = Pe + Pi ~ P e ■ 

Defining the momentum transfer by 

q = Pe - Pe = Pi - P' t 

and the average of the initial and final electron momenta by 

_ 1 , , 



(6.4) 
(6.5) 

(6.6) 



we can simplify Eq. (|6.3p to read 



01 



1 = 



dv e d»p e d v pi 



(2ith) v {2ith) v (2nh) 
1 



- T (2*H)8 — pi-q p-q- 



2m ? - 



P-Q /e(Pe)/i(Pi) 1 - /e(Pe) 



(6.7) 



Since T is a two-body scattering amplitude, its general 
form can depend upon both the square of the momentum 
transfer q 2 = q ■ q and the total center-of-mass energy 
W =p 2 /2m ei , where the relative momentum is given by 
p = m e j(v e — Vj), with m ei being the reduced electron- 
ion mass. It is the ^-dependence in T = T(W, q 2 ) that 
renders the integrals in Eq. (|6.7[) difficult to calculate 
because W depends explicitly on p.;. In Section 12 of 
Ref. [l[ , this calculation is performed to all orders in the 
Coulomb scattering. For the work here, we shall be less 
general and exploit the fact that the electron-ion mass 
ratio m e /m,i is very small (so the reduced mass m e i is 
almost equal to the electron mass m e ). We shall assume 
that the electron and ion temperatures are not orders 
of magnitude apart, a mild restriction in all practical 
applications, so that 

P e m e < fa rm . (6.8) 

Under these circumstances, the thermal average electron 
velocity is much larger than the ion velocity, and to a 
very good approximation |v e — Vj| = |v e |. Thus the 
quantum Coulomb parameter rj e i — e ei/47rh\~v e — Vj| that 
appears in the Boltzmann equation can be replaced by a 
Coulomb parameter that contains only the electron veloc- 
ity, rjei — > Zie 2 /4irH\v e \, where we have written ej = Zj,e. 



The size of this parameter is estimated by its thermal 
average, which we denote by an overline. We use the 
simple Maxwell-Boltzmann distribution to estimate this 
average. For this classical distribution, the thermal av- 
erage of 1/Vg is precisely m e /T e , and so 

where in the second equality we have used the previous 
definition (|3.4[) of the electron quantum Coulomb param- 
eter rj e together with the definition (|2.7p of the electron 
thermal wave length A e . Another way to write this is 

(6-10) 

where e H ~ 13.6 eV previously noted in Eq. (|1.11|) is the 
binding energy of the hydrogen atom. The result (|6.10|) 
demonstrates that rj e i is quite small for the elevated tem- 
perature range that concerns us. Hence the scattering 
amplitude in Eq. (|6.7p can be calculated in the Born ap- 
proximation [15l ]. 

T~T B (q 2 )=h^, (6.11) 
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a quantity that depends only upon the square of the mo- 
mentum transfer q 2 , and not on the center-of-mass en- 
ergy W. 

In the Born approximation, the initial ion momentum 
Pi appears only in the delta-function and phase-space 
density explicitly shown in Eq. (|6.7j) , and not in the am- 
plitude T B (q 2 ), and so the integration over this momen- 
tum variable can be carried out. If it were not for the 
delta-function factor, the p^-integration would simply en- 
tail 



(2-kK) 



-MPi) 



exp< -(3 j 



Pi 

2m,- 



(6.12) 



We shall often denote the magnitude of the momentum 
transfer by q = |q[, as we have done here. We now 
change the remaining two integration variables p' e and p e 
in Eq. (|6.7p to the variables p and q defined in Eq's. (|6.5p 
and (|6.6[) . a change that has a unit Jacobian. Since the 
electrons are described by the Fermi-Dirac distribution 
(12. 2|) , the Pauli blocking term in Eq. (|6.7|) can be written 
as 



Following the convention exhibited in Eq. (|2.5p , the 
species index i for ions implicitly includes spin degrees of 
freedom, and so the integration over a single /j(pi) pro- 
duces rii /$i . The delta- function in Eq. (|6.7p removes one 
of the components of the p^-integration, which is equiv- 
alent to supplying an extra factor of and retaining a 
Maxwcll-Boltzmann factor corresponding to the compo- 
nent of the momentum p; along the direction of q. Hence 



p»q p_q _r_ 

2m, 



in. 



in, 




Pi 



2m,i q 2 \ m 



nil _ q- 

2 , — p-q+- 



(6.13) 



l-/e(p + q/2) 



(6.14) 



g-/3ePe g xp< 



_Pe_ 

2m e 




/e(p + q/2) 



Using these results, and neglecting terms involving the 
very small ratios m e /mi and j3 e m e / j3 I irii 1 we find that 



dt 



= 2 E, 



d v p d v q 



(2itK) v (2nh) 
/ e (p - q/2) / e (p + q/2) exp 



r Pe 


- 2 , 1 2 


I 2m e 





Pi m, 

p • q exp^ 



p • q 



m e 
2m; 



i & I 



(6.15) 



where q = q/|q|, the variable p ± in the first exponent is 
the component of p orthogonal to the momentum trans- 
fer q, so that p = p ± + (p • q) q with p\ • q = 0. 



We can simplify the rate (|6.15[) by further exploiting 
the consequences of the very small ratio m e /?n,. In the 
order of magnitude estimates that follow, we will use the 
symbol (3 to designate the inverse temperatures of cither 
the electrons or the ions. This is possible because the 
temperature disparity is not very severe. We now see 
that the thermal distribution functions in the second line 
of Eq. (|6.15p restrict the size of the momenta to be of the 



order 



and 



^. (6.16) 



For the longitudinal component of the electron momen- 
tum, the form of the exponential in the last line of 
Eq. (|6.15p motivates the change of variables to 



P « = p,q+ 2^( 1 -A' ff 



(6.17) 
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Under this change of variables, the last line in Eq. (|6.15[) 
becomes 

m e ( (3 e \ 1 { ft m, _ /2 l 

(6.18) 

a term whose exponent restricts the size of the longitu- 
dinal component to be 




(6.19) 



This means that the second term in square brackets at 
the start of expression (|6 . 18[) . the term [m e /mi)q ~ 
yj mf//3rnf , is a factor y/m^Jrni smaller than the first 
term p!, . However, as we shall find, the first term inte- 
grates identically to zero, leaving the ostensibly smaller 
second term as the leading order contribution. To 
see this, we first note that the electron distributions 
/e(p T q/2) in Eq. (|6 . 1 5f> are functions of the dimen- 
sionless variables 

ft£ e (pTq/2) = A(p Tq / 2 )2. (6.20) 



Here, we must express the old variable p in terms of the 
new variable 



p' ee p ± +p' n q, , (6.21) 
or in terms of the vectors (16. 5p and (16. 6|) . 

Then from Eq. (|6. 22[) . we see that replacing the old vari- 
able p in Eq. (|6.20p by the new variable p ' incurs relative 
error of order 

K/m t )p' • q(lftp 2 ) ~ (m e /m,) 3/2 , (6.23) 

an error beyond the leading term that we retain. That is 
to say, we can simply replace 



/ e (p - q/2)/ e (p + q/2) -> / e (p' - q/2)/ e (p' + q/2) 



(6.24) 



This product is explicitly even in q, as are the remaining terms in the integrand, and consequently, the odd term 
in the prefactor of (|6.18[) integrates to zero. The energy rate (|6. 15p now reduces to 



dt 



/ e (p'-q/2)/ e (p' + q/2) exp 



2m P 



2m, 



q cxp 



Pi W^i - 1 2 

o p \\ 



(6.25) 



The integral over the momentum p' = p x + p n q contains v — 1 integrals from p^ and one integral from p!, . 

Now that the leading contribution has be extracted, we can make further reductions by omitting several more terms 
in m e /mi . In particular, we may now neglect the longitudinal part p!< = p' • q relative to q — |q| in the electron 

distribution functions / e (p' T q), which then become functions only of (p x =F q/2) 2 . In fact, since f> ± ■ q = 0, both 
electron distribution functions have the same argument, 



P e E e (p ± ±q/2) 



2m e 



1 , 
4*' 



(6.26) 



and their product becomes a simple square: / e (p ' — q/2) / e (p ' + q/2) = [f e (p± + q/2)] 2 - The longitudinal component 
p^ now appears only in the final factor of the integrand in Eq. (|6.25p . and we may therefore explicitly perform the 
integration over this part of the momentum, 



2ttH 



exp<{ 



.A. 

2m, 



ii 



(6.27) 
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where the ionic thermal wave-length is determined from Eq. (|2.7|) . We can now express the rate (|6.25[) as 



dt 



Pi 
-0 



- \ T B (q 2 ) \ 2 q 



Q l m l J {2nhy- 1 (2nh) 



Me [/ e (px + q/2) 



exp 



Pe 

2m P 



1 



Pa 



q 



(6.28) 



The integration over the momentum transfer q is damped at large values of q = |q|, because at such large values 

Pe > 1 



q — > oo 



/ e (Px+q/2) 



exp 



2m e 



Pa 



■ q 



g+/3eMe exp< ; 



f Pe 


- 2 i 1 2 




1 2m e 


PA + - A q 


} 



(6.29) 



Since the limit (|6.29[) constrains the integrand to small-g, this further supports the use of the Born Approximation 
(|6.1ip . which allows us to express the rate (|6.28[) as 



dt 



-f3 e e z m e h 2 u 2 [T e - T, 



Pa d^q 1 



{2nhy- 1 {2nh) l/ q 3 



f e (p ± + q/2) 



exp 



Pe 

2m e 



(6.30) 



where u> 2 is the sum over all the ionic species of the 
squared plasma frequencies, 



E ^ 2 = E 



ion species 



1 mi 



(6.31) 



Now that the clutter has abated, we can more easily 
study the nature of the parameters that enter into the 
v > 3 contribution of the energy exchange rate. Dividing 
the time derivative of the electron energy density by the 
the electron specific heat 3n e /2 gives the rate of the elec- 
tron temperature change already noted in (|5.7[) . namely 



dT 



(6.32) 



The integral of each momentum, with the normalizing 
denominator 2irh, gives a pure number times a factor of 
1/A e . Since n e ~ Aj w , and each factor of the momentum 
transfer |q| will produce a factor of H/X e , we conclude 
from Eq. (jOUf that 



7e_f ~ p e e 2 m e h 2 uj 2 Ag " (A e /fi) 
e 2 1 \ /fiw,' 



v Ae e / \ e 

where in the second line we have made use of 



A 



H 2 p e /m e and P e = 1/T e 



(6.33) 



(6.34) 



In a i/-dimensional space, the energy between two elec- 
trons a distance A e apart is, up to a constant, given by 



2 /K 



Hence the first factor in parenthesis in the last 



line above is dimensionless. Since Filo, is an energy, the 
second factor is also dimensionless. Thus the overall di- 
mension of 7 e j is that of the final factor w I} the correct 
dimension of an inverse time or rate. Although a factor 
of h appears here, it is canceled by the single factor of 
h that appears in l/Ag -2 in the v — » 3 limit, an so in 
this limit the rate is a classical quantity. However, as we 
shall see, the dimensional continuation method that we 
use produces logarithms, and a logarithm of h will ap- 
pear in the final result. Finally, we should note that the 
rate involves the first power of the ion density, a power 
that does not depend upon the spatial dimensionality v. 

The rate (|6.30[) for v > 3 diverges when v — > 3 + , a di- 
vergence that is canceled by the v — > 3~ limit of the rate 
for v < 3 that we compute in the next section. This latter 
rate involves purely classical dynamics. Thus it entails a 
wave-number k that comes from the Fourier transform of 
a potential which is the analog of the quantum momen- 
tum transfer q, but with q = ftk. With this replacement, 
the electron distributions would become / e (p' + ftk/2), 
but since only classical quantities appear in the forth- 
coming v < 3 contribution, in this part the electron dis- 
tributions must appear only as / e (p')- Thus, to separate 
out a part of the v > 3 Boltzmann expression for the rate 
that will combine in a simple fashion with the v < 3 con- 
tribution that we shall soon examine, we construct this 
part by making the replacement 



f e (p ± +q/2) 



exp 



Pe 2 

q 



8m 



/e(Px) 



exp 



Pe g2 



8m 



(6.35) 
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which exhibits the needed large q 2 damping given in the where 
limit (|6.29p . Accordingly, we decompose the rate of en- 
ergy transfer into a potentially singular part and a regular 
part, 

dS> dS> s d£>« 

^T^^T + ^T' (6 ' 36) 



dt 



-f3 e e z m e H 2 oj 2 ( T e - T, 



J e 2 



d u ~ 1 p ± d v q 1 
(2irft)"- 1 6XP 1 ~8^ q 



/e(pj 



exp 



2m, 



(6.37) 



and 



dt 



-l3 e e 2 m e h 2 u J 2 (T e -T^ / — 



/ e (p x +q/2) 



d Px d 3 q 1 e _ /3e/Je 



(27Tft) 2 (27Tft) 3 q 3 
1 



/e(Px) 



exp 



Pe (-2 1 2 

"2^7 \ Px ~A q 



(6.38) 



Here in the regular part, we have taken the v — > 3 limit since there is no impediment in so doing. 

The sing ular part (f6~37l) may be simplified by performing the g-integration. Passing to hyper-spherical coordinates 
gives 



d v q 
(2irh) u q 



1 

^3 ex P 



Pe 2 

q 



(2irhy Jo q q <x,, i s,i, 



& 2 



dx ( 8m 



0,-1 1 



(2-kK) v 2 J Q x \ 0, 



("-3)/2 



(2tt/i) 3 2^4 



2 x (3-f)/2 



(6.39) 



where f2j,_i is the area of a unit [y — l)-sphere embedded in a ^-dimensional space. In the second line above we have 
made an obvious change to a dimensionless integration variable x, and in the last line we have identified the resulting 
integral with a standard representation of the gamma function. We thus have 



dt 



Pe e 2 m e 

2h 



2 ^v — 1 

°' J2nf 

(2ith) v - 



ttA| 
4 



(3-1O/2 



Z/-3 



/.(Pi) 



exp<^ /?, 



2m e 



(6.40) 



As we shall see, the ^ < 3 contribution calculated in the next section contains the same integral over the transverse 
components p ± , so it will be convenient to perform this integral when we add these terms together. 

The regular part of the energy exchange rate may also be simplified since the integral over the electron distribution 
functions can be performed when v = 3. Namely, we pass to polar coordinates, with the angular integration simply 
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producing a factor of 2w, to obtain 

d 2 p ± 



/e(Pi) exp<^ (3, 



Hence, making the replacement 



for the first term, we have 



dt 



2ir(3 e H 2 



Pi 
2m 

2m P 



exp 



} 



exp{+/3 e [^-- Me ]J +1 



X 2 exp{-/3 e /x e } + 1 



j3 e He -» /3 e 



8m„ 



(3 e e 2 m e lo 2 1 



2fi 



dq 
<1 



exp 



8m c 



eX p{~&£} 

j + 1 exp{-/3 e /i e } + l 



/3 e e 2 m e uj 2 1 / \ - . 

— ^^A?C /; ^2 7 "' lnr 



exp{x - P e fl e } 



exp{— x) 



\ [exp-jz - P e H e } + l] 2 exp{-/3 e Me} + 1 J ' 

where the last line follows from a trivial change of integration variables and a partial integration. 



(6.41) 



(6.42) 



(6.43) 



VII. LENARD-BALESCU EQUATION: LONG-DISTANCE PHYSICS 

We turn now to calculate the leading order long-distance physics by working in spatial dimensions v < 3. This is 
done by employing the Lenard-Balescu equation with Fermi-Dirac statistics for the electrons and Maxwell-Boltzmann 
statistics for the heavy ions. For the spatially homogeneous system that we work with, the Lenard-Balescu equation 
with the appropriate Pauli-blocking reads 



dfe(Pe 

dt 



■—■T 



d v Vl d v k 



i J (2nh) v (2-k)" 



eei 



k 2 e(k, k • Vi 



k '^r /e(Pe) 1 _/,(p ' ] ~ /;(p;,k 



n5[ kv e - kv, 

dfe(Pe) 



(7.1) 



where the gradient d/dp e acts on everything to its right, and e(k, u) is the classical dielectric function for the plasma 
discussed in Appendix fXJ As was shown in Appendix C of Ref. [lj, the usual non-degenerate Lenard-Balescu equation 
is a formal limit of the Boltzmann equation. The same methods that were employed there may be used to show that 
Eq. (|7.ip is the corresponding long-distance equation when the electrons are degenerate and described by Fermi-Dirac 
statistics. If the electrons and ions are in equilibrium with themselves at temperatures T e and Tj, respectively, then 
their distribution functions f e and /, are given by Eq's. p.2| and (|2.ip . in which case the terms in curly braces can 
be written as 



V e -fiik- Vi) /i(pi) f e (Pe) exp<M, 



Pe 

2m P 



(7.2) 
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Because the delta- function equates k • v e with k • Vj, the factor (|7.2p and with it the right-hand-side of Eq. (|7.ip 
vanishes when the electrons and ions are in thermal equilibrium with a common temperature T, the electrons being 
described by a Fermi-Dirac distribution and the ions by a Maxwell-Boltzmann distribution. This confirms the validity 
of Eq. (|7.ip . the Lenard-Balescu equation with Pauli-blocking for degenerate electrons. 

Using (|7.2p with the ions at the same inverse temperature /3 7 = 1/Tj, and upon integrating the total derivative 
d/dpe by parts, we can express the energy exchange rate as 



dt 



2 £, 



d v p e (Tpi d v k 



k 2 e(k, k • Vi) 



7T <5^k-v e 



k- v, 



(fie k • V e - (3j k • Vj) fi(pi) /e(p e ) exp j 



2m, 



(7.3) 



We have placed a "less-than" superscript on the left-hand 
side of (|7.3p to remind ourselves that the calculation is 
performed in v < 3 using the Lenard-Balescu equation. 
The distribution functions constrain the velocities of the 
ions and the electrons to be of the order Vi ~ yjTj/rrii 
and i; e ~ y/T e /m e , respectively. Since an ion mass is 
so much greater than that of an electron, and the tem- 
perature disparity is never excessively large for cases of 
interest, the ions move much slower than the electrons, 

Vi <C v e , (7.4) 

This restriction also follows from the previous condition 
(|6.8p . the condition that (3 e m e <C f3 I rrii. To compute the 
rate (|7.3p . we first decompose the electron momentum 
into perpendicular and longitudinal components relative 
to the direction specified by k, so that p e = + j?||k 
with k • p^ =0 and = k • p e = m e k • v e . The delta 



function in Eq. (|7.3p can be used to remove the paral- 
lel component p\\ of the electron momenta integration. 
Since dp\\/(2irh) — (m e /2nhk) rf(k • v e ), the use of this 
delta function produces a factor (m e /2Trhk) and makes 
the replacement p\\ — ► m e k- Vj. In view of the limit (|7.4p . 
this replacement makes p\\ much smaller than the mag- 
nitude of the perpendicular components of the electron 
momenta p ± . Hence, we can simply replace p — > p ± in 
the remainder of the integrand. We shall also find it con- 
venient (the heart hath its reasons), to insert a factor of 
unity in the form 



1 = / dv 5(u-k-Vi) , (7.5) 



which allows us to express Eq. (|7.3p as 



dt 



= -2/3, 



2nh 

d v k , 
k 



d"~ 



{2nhy- 1 

„,2 



dv 



/e(p±)J exp<{ (3, 
d v p. 



2m„ 



k 2 e(k, vk)\ 



2 



Pxet 



{2-Khy 



fi(Pi)s( 



k- 



(7.6) 



This integral can be further simplified by taking advantage of the analytic properties of the dielectric function 
e(k,w), discussed in some detail in Appendix lAl Repeating Eq. (|A16|) here for convenience, we see that a considerable 
portion of the integral simplifies because 



\k 2 e(k lV k)\ 



2 Z-^ 



{2-Khy 



,f l {p i )8{ii-^i 



i 

~2i 



k 2 + Kl + Fj(v) k 2 + K 2 e +F I {- 



(7.7) 



This result is the (unemotional) reason that the factor of unity in the form displayed in Eq. (|7.5|) was inserted in 
the integrand. Here n e is the electronic contribution to the Debye wave number, including the effects of Fermi-Dirac 
statistics, as expressed by Eq. (IA6|) . while Fj is a complex- valued function defined by Eq. (|A9|) . It is important to 
realize that the simplification (|7.7p only occurs when the ion species are summed over. The function F^z) is analytic 
over the upper half of the complex z-plane, and has the asymptotic behavior 
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where w ; is the total ionic plasma frequency defined above in Eq. (|6 . 3 1 1) . Since an explicit odd factor of v appears in 
the integrand, we can write the resulting integral over v in Eq's. (|7.6[) in the form 



+ OC 



dv — 

2i 



1 



1 



k 2 + k 2 + Fj (v) k 2 + k 2 + Fj(-v) 



= lim % 



+v 



dv 



-V 



k 2 + k 2 + F(v) 



(7.9) 



The delta-function in Eq. (|7.7[) removes the longitudinal components of the ionic momenta, leaving a Maxwell- 
Boltzmann factor involving the velocity v. Hence the left-hand-side of the integrand in Eq. (|7.9p is damped in a 
Gaussian fashion for large |u|. This rapid damping results from a cancellation between the terms with F,(v) and 
Fj^—v) that only happens when the two velocities are exactly the negative of one another. Hence, when we simplify 
the integrand by taking advantage of the odd prefactor as was done on the right-hand side of Eq. (|7.9[) . we must 
integrate between the exact same negative and positive limits, between — V and +V , and only afterward take the 
limit V — > oo. Since F I (z) is analytic in the upper-half z-plane, the integral (|7.9p may be evaluated by contour integral 
techniques. Let C v be a semicircle of radius V centered at the origin of the complex z-plane, with an orientation that 
starts at +V and ends at —V. We can traverse a closed circuit by moving from — V to +V along the real axis, with 
the circuit completed back to —V by traversing C v . The contour integral around this closed circuit vanishes since it 
contains no interior singularities, 







dz 



k 2 + K 2 +F I { Z ) 



dv ■ 



-v 



dz- 



k 2 + K 2 e + Fj (v) J Cv k 2 + K 2 e + Fj (z) 



(7.10) 



Hence the integral (|7.9p is equal to the negative of the integral over the semicircle C v starting at +V and ending at 
— V. We can now take the limit V — * oo and use the asymptotic form (|7.8p for F, along C v . Since 



dz z — i V 2 



d9e 2lB 



= 0, 



(7.11) 



the leading term of the expansion of the denominator in Eq. (|7.9p yields a vanishing contribution, and the only 
non-vanishing term is given by first-order term with Fj(z) replaced with its asymptotic form (|7.8p : 



lim i 



+V 



dv 



1 



k 2 + Fj(v) 



[k 2 + «2 

7T 



lim 

V^oo 







/ dz z 

ICv 


z 2 . 



[k 2 + KlY 

Upon passing to hyper-spherical coordinates to perform the ^-integration, we now arrive at 



(7.12) 



dt 



Pe e 2 m e 



1V-X 



p± 



2 ^v — 1 



/e(pj 



exp<^ (3, 



PI 



(27T)" 



■ v-X 



dk k ■ 



2m e 
1 



[k 2 



(7.13) 



K7. 



Changing variables by k = t x l 2 K e places the fc-integration in the form of a standard representation of the Euler Beta 
function [16[, and so we have 



dk- 



k" 



[k 2 



1 r f^±i 

2 e 



r(2) 



(7.14) 



Finally, we are now able to express the v < 3 form of the electron-ion energy exchange as 



dt 



(3 e e 2 m e 
2ft 



2 "j/ — 1 / K e V 

1 T2W \2^J 



(27T) 



1 



(271-ft) 1 - 



/.(Pl) 



exp-^ (3, 



2m P 



(7.15) 
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VIII. ADDING THE RATES 



The sum of the singular part (|6.40[) for the v > 3 contribution to the electron-ion energy exchange rate and the 
v < 3 part (|7.15|) that we have just computed is 



dt dt 



& 777- g n ^~^>i/— 



2H 1 (2tt) 3 



?^( T -- r ')yfe[AM"<H* 



(27rft) 



2m„ 



.o \ (3-v)/2 f / o\ / 2\2\( I/ - 3 )/ 2 / ,\ / o 



16tt 



(8.1) 



As must be the case, the expression in the final curly braces above is finite in the v — > 3 limit. To extract this limit, 
we use 



v -> 3 



Z/-3 



J/-3 



7 , r 



3-i/ 



3-i/ 



7 , r 



z/+ 1 



l-(l- 7 ) 



3-f 



to evaluate the v — > 3 limit of the last line in Eq. (|8.1 



i/-3 



7 



„2\2\ (^-3)/2 
«„ A, 



16tt 



3-1/ 



- 1 



hi 



16tt 



-7-1 



(8.2) 



(8.3) 



Since this last factor is finite in the v — > 3 limit, we may now take the v — > 3 limit of all the other quantities in 
Eq. (|8.ip . The integral over the perpendicular momenta in this limit was evaluated previously in Eq. (|6.4ip . and so 
we now have 



d£> s dS< P e e 2 m e ujj 1 1 



In 



16tt 



-7-1 



^ 9t 2h 2tt 2 A2 exp{-/3 eA i e } + 1 

To this we must add the remaining finite part (|6.43[) of the v > 3 contribution, namely 



T e - T, 



d£> R P e e 2 m e u 2 1 



dt 2h 2 

Recalling that we have defined [Eq. (|5.1[) ] 



exp{x - /3 e [i e } exp{-a;} 



[exp{a; - p e fi e } + l] 2 exp{-/3 e fi e } + 1 



(8.4) 



(8.5) 



d£ eI 
dt 



■C, 



■6) 



we have now calculated the rate coefficient to leading order in the plasma coupling and to all orders in the electron 
fugacity z e = e^ eAte , 



/3 e e 2 m e lu 2 1 



1 1 



2h ir 2 \ 2 e I exp{-/3 e fi e } + l 2 



In 



16tt 



7-1 



+— / In a: 



exp{x - /3 e /i e } exp{-x} 



[exp{x — /3 e Af e } + l] 2 cxp{-/3 e + 1 
By expanding the denominators, it is easy to check that 

exp{x - p e fi e } exp{-ir} 



•7) 



dx In. x 

[exp{x - p e fi e } + l] 2 exp{-/3 e fi e } + 1 
which is an expansion in powers of the electron fugacity z e = e 13 "^" 



]T(-1)' +1 ln{? + l}e(' +1 



(8.8) 
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To place this result in a form that is easily compared to that of BPS [l[, we use the definition (|2.7[) of the thermal 
wavelength and a slight manipulation to write 



C eI 



2ttV 2ir I A? 6 



exp{/3 e (i e } + 12 



In 



16tt 



7-1 



if>l)< +1 ln{/ + l}e^l 



(8.9) 



In the dilute limit in which Maxwell-Boltzmann statistics apply, the fugacity exp{/3 e /it e } is very small. The number 
density approximation (|2.12|) gives 



A 3 . 



1+ 2^ e 



PeUe 



and we see that keeping the first correction in the fugacity yields 

1 



ujf / I3 e m e 



2vr V 2?r 



(/3 e e n e ) 



1-1- 



2 3 / 2 





1 




' 16tt 1 






ln<j 


-7-1 




2 







+ -e^ In 2 



Again remembering the fugacity approximation (|2.17p . which we repeat here, 

1 



(3 e e n e 



1 



2 3 /2 



and the definitions 



we find that 



a: 



2 _ 2?r h 2 (3 e 



and 



e n e 



16tt 8T„ 2 



1+ 2W2 e 



and thus 



r ~ — 



/ (i e m e 



2tt V 2?r 



(P e e 2 n e 



1-1 



2 3/S 



hi 



87? 



(8.10) 



(8.11) 



.12) 



.13) 



.14) 



- In 2 + 

2 2 5 / 2 



.15) 



We may use e^"^' ~ X^,n e /2 inside the curly braces of 
(|8.15p . It is easy to confirm that with the neglect of the 
fugacity corrections, this is in agreement with Eq. (12.12) 
of BPS [l[ after that equation is corrected as mentioned 
in the Introduction. 
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APPENDIX A: THE DIELECTRIC FUNCTION 



In Section IVII1 the calculation of the rate in v < 3 using the Lenard-Balcscu equation made extensive use of the 
plasma dielectric function and its various properties. The classical dielectric function for a collisionless plasma is 
discussed in Ref. [TtJ for example, and the form of the result that we shall use reads 

e(k, U ) = 1 + V f| / *^ J— - k • A /b(pfc ), (Al) 
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with the prescription 77 — > + defining the correct retarded response. The degenerate electrons are described by the 
thermal Fermi-Dirac distribution (|2.2p . so 

Q e f3 c (E e -ti c ) 

= -/3 e k-v e / e (p e )[l-/ e (p e )] . (A2) 



On the other hand, the ions are described by the Maxwell-Boltzmann distribution (|2.ip . which is simply the large 
chemical potential limit — f3fi 3> 1 of the Fermi-Dirac distribution. In this limit the Pauli blocking term is removed, 
[1 - /(P)] -> 1) and so 

k~/ i (p i ) = -Ak-v i / i ( Pi ). (A3) 



For the real plasma considered in the text, the ions equilibrate to a common temperature T 1 = l//3 7 ; however, for 
the purposes of this Appendix, we shall take each ion species i to have an individual inverse temperature For 
degenerate electrons and Maxwell-Boltzmann ions, the dielectric function (|A1|) may therefore be expressed as 

<k,v) = 1-^-2 f ^^ /e (p e )[ We ( Pe )] 

k J (2Kh) v uj — k-v e + i?/ 



AC ' 1 d ""> k ' Vl /.(p.). (A4) 



k 2 J (2irh) u uj - k-Vi + ir\ 



The factor of two in the electron contribution arises from 
a sum over the two spin components of the electron. 

The dielectric function in the Lenard-Balescu equation 
has the functional form e(k, v ■ k), with the speed |v| 
much less than the electron thermal velocity. Hence, in 
the electron contribution to the dielectric function, the 
magnitude of w = k • v is much less than the typical 
magnitude of k • v e , and we can use the uj — > limit in 
which 



Pee 1 



■ 2 



d?p e 



k- 



(2irhy uj - kv e 



-/e(Pe) [1 - /e(Pe)] 



k 2 



(A5) 



where 



/ e (p e ) [1 - / e ( Pe )] (A6) 



defines the electron contribution to the squared Debye 
wave number, including the effects of Fermi-Dirac statis- 
tics which are explicitly exhibited by the Pauli blocking 
factor [1 — / e (Pe)]- From the form (|2.2p of the thermal 
Fermi-Dirac distribution / e (p e ) and the definition (|2.5p 
of the number density, we see that 



= e 2 /3 e 



(A7) 



Remembering the structure of the grand canonical en- 
semble, the derivative that appears here is the thermal 



average of the fluctuations about the mean particle num- 
ber. In the limit of Maxwell-Boltzmann statistics, the 
derivative simply reproduces the particle number density, 
corresponding to the fact that classical statistics have a 
Poisson distribution. Multiplying Eq. (|A5|) by k 2 and us- 
ing Eq. (|A5j> for the electron contribution, we can write 



(A8) 



k 2 e{k,uj) = k 2 + K 2 e +F I {uj/k) 7 
where we have defined the function 



k- 



k-Vi + IT] 



/i(Pt 



(A9) 



This is almost same function F defined in Ref. ex- 
cept that here we have handled the electron contribution 
separately. 

Superficially it appears that F, contains wave-vector 
dependence through the terms k • v, in the integrand 
of (|A9j) ; however, since we are integrating over all val- 
ues of Vj, the wave- vector direction k cancels in Fj. As 
our notation suggests, F T (v) is indeed only a function of 
v = |v|. Furthermore, because of the i^-term with r\ > 
in the denominator, the function F t is analytic in the 
upper complex u-plane. 

In evaluating the integral (17. 12[) in the text, we require 
the large- v behavior of (|A9|) . Since the numerator of 
the integrand in Eq. (|A9[) is odd, we can expand the 
denominator to find the leading u-behavior 
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where 

* — '* rrii 

is the sum of the squared ion plasma frequencies. To obtain this result, the Gaussian integral in Eq. (|A10|) may be 
calculated directly, or more elegantly, one may use Eq. (|A3[) to replace (k • v-j)/, with a derivative of ft, after which a 
partial integration yields Eq. (|A10j) . 

Finally, we derive a dispersion relation that will be quite useful in evaluating Eq. (|7.6p in Section IVIII Applying 
the relation 

Im - - — = -n6(x) (A12) 
for 77 — > + in Eq. (|A9[) allows us to express the imaginary part of Fj in the form 

Imi^)=E.fte* J -^-fifajvnSiv-k-Vi). (A13) 

From this, we can find the imaginary part of the inverse of the dielectric function: 

1 Imfc 2 e(/c,v-k) _ Im Fj(v) 

m fc 2 e(fc,vk) _ ~ |fc 2 e(A:,v-k)| 2 ~ ~ \k 2 e(k, v k)| 2 

- - \k 2 <k \ -w e <Jw¥ J ' (p ' ) " S{v ~ t - T|) ■ (A14) 

Since the numerator in the integrand (|A9p is odd, under complex conjugation we have 

F(-v) = F{+v)* . (A15) 

Hence, using Eq's. (|A8|) and (|A15|) can write 



|jfc 2 e(fc, 

1 r 1 1 



2i {k 2 + K 2 e + F(v) k 2 + K 2 e + F(-v) 
I 



(A16) 



APPENDIX B: SCATTERING CORRECTIONS 

The electron-ion energy exchange rate computed in 
Section 12 of BPS l| was performed under quite gen- 
eral conditions, with no restriction on the masses, num- 
ber densities, or temperatures of the plasma components, 
except that the plasma be fully ionized, non-degenerate, 
and weakly to moderately coupled (all mild restrictions 
in a hot, low-Z plasma). Because of its generality, this re- 
sult, which we shall present momentarily, is rather com- 
plicated. However, for most practical calculations, we 
can work in the high-temperature extreme-quantum limit 
and take advantage of the small electron-to-ion mass ra- 
tio. Under these conditions, we can use the Born approx- 



imation for the two-body scattering amplitude, and the 
rate coefficient collapses to the simple expression, 

u 2 lp e m e 2 1 I" f 8T 2 1 " 

c - = 27v^r (/3eene) 2rte/ 7 ' 

(Bl) 

as previously quoted in (|1.5[) . The purpose of this Ap- 
pendix is to find the subleading quantum correction to 
Eq. (Bp . 

As noted in BPS, this subleading correction is of order 

f&~27eV/T e , (B2) 

which, for most applications that we have in mind, is 
quite small. However, the leading electron degeneracy 
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effects are of order z e ~ n e a 3 , (27 eV/T e ) 3 / 2 , which can 
be comparable in size to the subleading quantum correc- 
tion. Both corrections are small compared to the leading- 
order contribution (|B1[) . Therefore, in this Appendix, 
we can work in the non-degenerate limit, since degen- 
eracy effects on top of the subleading quantum effects 
are smaller still. Consequently, our starting point will 
be the non-degenerate, but otherwise rather general, ex- 
pression for the rate derived in Section 12 of BPS. We 
shall simplify this rate in favor of the more realistic case 
of light electrons and heavy ions, exhibiting this result in 
Eq. (|B24[) . an expression that is valid to all orders in the 
quantum parameter r\ e i . Although this expression is clear 
and compact, for the purposes of this paper, however, we 
only need the subleading rj^ term. This subleading cor- 
rection is displayed in Eq. (IB33[) . 

The strength of the quantum effects associated with 
the scattering of two plasma species a and b is character- 
ized by the dimensionless parameter 



Vab 



e a e b 



(B3) 



The extreme quantum limit, where formally ft — > oo, is 
given by fj a b 0; while the extreme classical limit, where 
formally ^. — > 0, is given by fj a b — > oo. The former case is 
equivalent to the Born approximation. In Section 12 of 
BPS, the energy exchange rate from an arbitrary plasma 
species a to another species 6, 



d£ a b 
dt 



■ C a b (T a — Tb) 



(B5) 



was computed to all orders in the two-body quantum- 
scattering parameter fj a b. It was found that the rate co- 
efficient can be written as a sum of three terms, which, 
in the notation of BPS, reads 



Cab — C<b. n + [Cab,s + C a 



7AQ 



(B6) 



where the square of the thermal velocity in this expres- 
sion is defined by 



v\ = 

v ab 



l 



l 



f3 a m a (3 b m b 



(B4) 



where the last two terms have been grouped together 
for later convenience. These three terms are given by 
Eq's. (12.31), (12.25), and (12.50) respectively in BPS: 



ib,R 



K a K b 

2tt 



(3gm a 
2tt 



1/2 



p b m b \ 
2tt J 



dv v 2 e 



^{Pa.m a +f3 b m b )v 2 J_ 



F{v) 



^TT ptotal(^) 



hi 



F(v) 
K 2 



(B7) 



ib,s 



2 (I3am a (3bmb) 1/2 



(Pam a + Pbm b ) 



1 



3 /2 \2ir 



\i/2 



hi 



K 



e a e b 



4tt Am ab V 2 b 



2 7 



(B8) 



and 



'-ab 



1 2 2 ((3 a m a Pbm b ) 1/2 



{(3 a m a +Pbmb) 3 / 2 \27r/ J 



3/2 



d(e 



-C/2 



-111 



Vab 
(1/2 



(B9) 



Since we are using Maxwell-Boltzmann statistics 
throughout this section, the Debye wave number of any 
plasma species, including electrons, is here determined 
by n 2 = (3be 2 nb- The complex- valued function F(v) is 
defined by 



F(v) = / 



du 



Ptotal(w) 



II] 



where ptotai(v) is the spectral weight, 



Ptotal(w) = ^2 Pb(v) 



(BIO) 



(Bll) 



with 



Pb(v) 



(i b m b 
~2V 



■ v exp 



--fymbv' 



(B12) 



This is similar to the function Fj introduced in the pre- 
vious Appendix, except here the sum extends over all 
plasma species, including electrons. Note the depen- 
dence in Eq's. (|B7|) and (|B8[) on an unspecified parame- 
ter K, with the only restriction being that K has units 
of a wave-number. This is an artifact of the calcula- 
tional procedure, and it was shown in BPS that the total 
rate C ab is indeed independence K, as this parameter 
cancels in the sum between Eq's. (|B7[) and (|B8[) . As 
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a matter of technical convenience, we will henceforth set 
K = K e throughout the rest of this Appendix [the simpli- 
fied result (|B18j) only holds under this condition] . Finally, 
we should note that the reduced mass is determined by 
l/m Q b = l/m a + l/mb, and 



1 dT(z) 



r(z) dz 



(B13) 



is the logarithmic derivative of the gamma function. 



Specializing to electrons and ions (in which a = e and 
b = i), we can employ Eq. (|B12j) to write Eq. (|B7|) in the 
form 



i/o \ 1/2 



2tt 



(B14) 



This expression greatly simplifies since the electron is so much lighter than the ions. In virtually all practical ap- 
plications, the electron and ion temperatures are never excessively disparate, and we can therefore impose the mild 
restriction 



/3 e m e < fa-mi 



(B15) 



We will refer to the condition (|B15[) as the m e — > limit, and the ratio f3 e m e / piffii can then be used as a small 
dimensionless expansion parameter. For example, to leading order in this parameter we find 



Ptotal(f) 



= 1+0 



Am* 



1/2 



which allows us to express Eq. (|B14[) as 



(3 e 1Tl c 

2vr V 2?r 



1/2 



cZd w i^fu) ln 



F(v) 



l + O 



Pern, 
f3 i m l 



1/2 



(B16) 



(B17) 



We have omitted the exponential in the integrand of (|B14|) . since the function F{v) provides enough convergence at 
large values of v to allow the m e — ► limit to be brought inside the integral. The analytic properties of F(v) allow 
us to perform the v integral in Eq. (|B17|) using contour integral techniques, in much the same manner as we did in 
the discussion following Eq. (|7.9p in the text. The result is Eq. (12.44) of BPS, which reads 



f3 e m e <C fami 



14 
2 2tt 



f3 e m t 
2tt 



1/2 



u>7 



(B18) 



a much simpler expression indeed. For the electrons and ions we are considering, we can also drop the term of order 
/3 e m e / Piirii in Eq. (|B3[) . allowing us to express the quantum parameter of Eq. (|B4[) as 



Ve 



1 — = Zl ^ 



On occasion, we will express the quantum parameter in terms of the binding energy of the hydrogen atom 

€h= 2 



e 

47T 



h 2 



13.606 eV 



so that 



-2 ry2 

Vei = Z i 



26, 



(B19) 



(B20) 



(B21) 



We see that r] ei <C 1 when T e reaches the keV scale, illustrating that quantum corrections are important at high 
temperatures. Finally, we can drop terms of order f3 e m e / fiirrii in the leading coefficients of Eq's. (|B8[) and (]F39|) . 
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thereby allowing us to write 

j3 e m e < Pirrii 



nC , riAQ _ 
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Unlike Eq. (|B18|) . which only holds for the sum over ion components, the result (|B22|) actually holds component by 
component. Performing the £ integral for the last term in Eq. (|B22|) gives 
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The rate C eI is given by adding Eq's. (|B18[) and (|B23|) . which can be written as 

j3 e m e < Pimi : 



Ce 



Kg / (3 e m ( 



2vr V 2tt 



1/2 
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(B25) 



This expression is accurate to leading and next-to-leading order in the plasma coupling, and to all orders in fj e i, with 
no restriction on the temperature (apart from requiring the mild constraint (|B15j) and that the plasma coupling be 
small). 

With the aid of Eq. (10.17) of BPS, 
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^(l + zryC 1 / 2 ) +7 = E 



1 ff 



k=l 



k k 2 Q + rj 2 



and writing 



a partial integration now gives 



C + fj 2 /k 2 dC, 
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(B27) 



(B28) 



In the latter form, we can easily extract the leading order term r] e i, since we can use the limit 
fjei -* : 



dCe~ c/2 lni( 



k 2 



rf(C/2)e- c/2 (ln{C/2} + ln2 
7 + In 2 . 
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Using now 



OO 



fc=l 



(B30) 



27 



and 



C'(3) = -^T p lnfc = -0.19812- ■■ 



fe=i 



we can express the subleading quantum correction as 

Mfjei) ~ C |c(3) ln{^}-7 -2C'(3)}. 

Using Eq. (|B21ll . we can express the rate (|B24|) to leading and next-to-leading order in fj e i as 
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The correction is of order e H /T e , which can be of comparable size to the leading degeneracy correction. 



[1] L. S. Brown, D. L. Preston, and R. L. Singleton Jr., Phys. 
Rep. 410, 237-333 (2005); arXive: |physics/050i084| 

[2] For a more detailed treatment of the method of dimen- 
sional continuation see: R. L. Singleton Jr., LA-UR-06- 
6738; arXive:0706.2680. 

[3] L. S. Brown, Ph ys. Rev. D 62, 045026 (2000); arXive: 
lphysics/9911056^ 2. 

[4] H. Brysk, Plasma Phys. 16, 927 (1974). 

[5] See, for example, H. Brysk, P. M. Campbell, and P. Ham- 
merling, Plasma Phys. 17, 473 (1975). 

[6] See, for example, C-K. Li and R. D. Petrasso, Phys. Rev. 
Lett. 70, 3059 (1993), where a similar minimum distance 
is used for the motion of a charged particle of given speed 
rather than a thermal ensemble of particles that we con- 
sider. 

[7] This is a consequence of Carlson's theorem, Upsala thesis, 
1914. The theorem is proven, for example, in E. C. Titch- 
marsh, The Theory of Functions, 2nd Ed., Oxford Univ. 
Press (London, 1960), See's. 5.8 and 5.81. 

[8] L. S. Brown, Quantum Field Theory, Cambridge Univ. 
Press (Cambridge, 1992). 



[9] J. Hubbard, Proc. Roy. Soc. A 261, 271 (1961). 
[10] E. A. Frieman and D. C. Book, Phys. Fluids 6, 1700 
(1963). 

[11] H. A. Gould and H. E. DeWitt, Phys. Rev. 155, 68 
(1967). 

[12] O. Aono, Phys. Fluids 11, 341 (1968). 

[13] R. L. Liboff, Kinetic Theory 2nd. Ed., John Wiley and 

Sons, Inc. (New York,1998). 
[14] Liboff uses the notation 7o for the Boltzmann kernel, I for 

the Lenard-Balescu kernel, and Io for the renormalization 

term. 

[15] The validity of using the Born approximation for the scat- 
tering amplitude is discussed in more detail in Ref. [J], 
particularly in the work leading to its Eq. (12.54). 

[16] See for example, § 6.2.1 of M. Abramowitz and LA. Ste- 
gun, Handbook of Mathematical Functions, Dover Publi- 
cations, New York, 1972. 

[17] E. M. Lifshitz and L. P. Pitaevskii, Physical Kinetics, 
Volume 10 of the Course of Theoretical Physics, § 29, 
Pergamon Press, Oxford, 1981. 



